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Abstract. The dynamic algorithm to compute a Grobner basis is nearly 
twenty years old, yet it seems to have arrived stillborn; aside from two initial 
publications, there have been no published followups. One reason for this may 
be that, at first glance, the added overhead seems to outweigh the benefit; 
the algorithm must solve many linear programs with many linear constraints. 
This paper describes two methods that reduce both the size and number of 
these linear programs. 

1. Introduction 

Since the first algorithm to compute Grobner bases was described by [8], they 
have become a standard tool for applied, computational, and theoretical algebra. 
Their power and promise has stimulated a half-century of research into computing 
them efficiently Important advances have resulted from reducing the number of 
pairs considered |5J [TSl [TU] , improving the reduction algorithm |31l [HI [5] , 
and forbidding some reductions [321 [TJ [T3] . 

The Grobner basis property depends on the choice of term ordering: a basis can 
be Grobner with respect to one term ordering, but not to another. Researchers 
have studied ways to find an ordering that efficiently produces a basis "suitable" for 
a particular problem [29, 30j, and to convert a basis that is Grobner with respect 
to one ordering to a basis that is Grobner with respect to another |161 . 

Another approach would be to begin without any ordering, but to compute both 
a basis and an ordering for which that basis is Grobner. Such a "dynamic" algorithm 
would change its ordering the moment it detected a "more efficient" path towards 
a Grobner basis, and would hopefully conclude with a smaller basis more quickly. 

Indeed, this question was posed nearly twenty years ago, and was studied both 
theoretically and practically in two separate papers |20l [9]. The first considered 
primarily questions of discrete geometry; the dynamic algorithm spills out as a 
nice application of the ideas, but the authors did not seriously consider an imple- 
mentation. The second described a study implementation that refines the ordering 
using techniques from linear programming, and focused on the related algebraic 
questions. 

Aside from one preprint [IS], there has been no continuation of this effort. There 
certainly are avenues for study; for example, this observation at the conclusion of 

m- 

In a number of cases, after some point is reached in the refining of 
the current order, further refining is useless, even damaging. . . . The 



2000 Mathematics Subject Classification. 13P10 and 68W30. 



DYNAMIC GROBNER BASIS ALGORITHMS 



2 



exact determination of this point is not easy, and an algorithm for 
its determination is not given. 

An example of the damage that can occur is that the number and size of the linear 
programs grow too large. The temporary solution of [5] was to switch the refiner 
off at a predetermined point. Aside from the obvious drawback that this prevents 
useful refinement after this point, it also forces unnecessary refinement before it! 
This can be especially costly when working with systems rich in dense polynomials. 

This paper presents two new criteria that signal the refiner not only to switch 
off when it is clearly not needed, but also to switch back on when there is a high 
probability of useful refinement. The criteria are based on simple geometric insights 
related to linear programming. The practical consequence is that these techniques 
reduce both the number and the size of the associated linear programs by significant 
proportions. 



2. Background 

This section lays the groundwork for what is to follow, in terms of both termi- 
nology and notation. We have tried to follow the vocabulary and notation of jS], 
with some modifications. 

Section [2~T] reviews the traditional theory of Grobner bases, inasmuch as it per- 



tains to the (static) Buchberger algorithm. Section 2.2 describes the motivation and 
background of the dynamic algorithm, while Section |2.3| reviews Caboara's spec- 
ification. Section |2.4| deals with some geometric considerations which will prove 
useful later. 



2.1. Grobner bases and the static Buchberger algorithm. Let m,n € N + , 

K a field, and R = F [x%, . . . , x n ]. We typically denote polynomials by /, g, h, p, 
q, and r, and the ideal of R generated by any F C R as (F). Following [9 , we 
call a product of powers of the variables of R a term, and a product of a term 
and an element of K a monomial. We typically denote constants by letters at the 
beginning of the alphabet, and terms by t, u, v. We denote the exponent vector of 
a term by its name in boldface; so, if n — 4 and t = x\ x^xf*, then t = (2, 0, 1, 20). 

An ordering a on the set T™ of all terms of R is admissible if it is a well- 
ordering that is compatible with divisibility and multiplication; by "compatible 
with divisibility," we mean that t \ u and t ^ u implies that t < a u, and by 
"compatible with multiplication," we mean that t < a u implies that tv < a uv. We 
consider only admissible orderings, so henceforth we omit the qualification. 

We write T for the set of all term orderings, and denote orderings by Greek 
letters fj,, a, and r. For any p 6 R we write lt CT (p) and lc CT (p) for the leading term 
and leading coefficient of p with respect with a. If the value of a is clear from 
context or does not matter, we simply write It (p) and lc (p). For any F C R, we 
write It, (F) = {lt CT (/):/€ F}. 

Let / be an ideal of R, and G C I. If for every pel there exists g g G such that 
It (g) | It (p), then we say that G is a Grobner basis of I. This property depends 
on the ordering; if a, t G T, it is quite possible for G to be a Grobner basis with 
respect to a, but not with respect to r. 

It is well known that every polynomial ideal has a finite Grobner basis, regardless 
of the choice of ordering. Actually computing a Grobner basis requires a few more 
concepts. Let f,p,r E R. We say that p reduces to r modulo /, and write 
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inputs: 

• F C R 

• a G T 

outputs: G C i?, a Grobner basis of (F) with respect to cr 

do: 

(1) Let G = {}, P = {(/,0) :/£F} 

(2) while P f 

(a) Select (p, g) G P and remove it 

(b) Let r be a remainder of spoly (p, q) modulo G 

(c) if r^O 

(i) Add (g, r) to P for each g G G 

(ii) Add r to G 

(3) return G 

FIGURE 2.1. Static Buchberger Algorithm 



p — 1-> r, if there exist a E K and i G T" such that p — a£/ = r and It (r) < It (p). 
Similarly, we say that p reduces to r modulo G, and write 

V — > r , 
G 

if there exist . . . , ig} C {1, . . . , ^G} such that p — > r 1; rx — ^ r 2 , . . . , rg_i — >• 

= r. If there no longer exists g G G such that It (5) divides a It (r), we call r a 
remainder of p modulo G. 

Proposition 1 (Buchberger's Characterization, G is a Grobner basis of I if 
and only if the S-polynomial of every f,g G I\ {0}, or 

sP oiy </,») - ic (9 ) . lc -" t 1 ;{} s '"»» . / - ic </) . -■" t 1 tf» i "fa>» . 9 , 

reduces to zero modulo G. 

For convenience, we extend the definition of an S'-polynomial to allow for 0: 

Definition 2. Let p G R. The ^-polynomial of p and is p. 

Buchberger's Characterization of a Grobner basis leads naturally to the classical, 
static Buchberger algorithm to compute a Grobner basis; see Algorithm |2.1| 
which terminates on account of the Hilbert Basis Theorem (applied to (It (G))). 
There are a number of ambiguities in this algorithm: the strategy for selecting 
pairs (p, q) G P, for instance, or how precisely to reduce the ^-polynomials. These 
questions have been considered elsewhere, and the interested reader can consult the 
references cited in the introduction. 



Remark 3. Algorithm 2.1 deviates from the usual presentation of Buchberger's 
algorithm by considering .^-polynomials of the inputs with 0, rather than with each 
other. This approach accommodates the common optimization of interreducing the 
inputs. 
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2.2. The dynamic algorithm. Every admissible ordering can be described using 
a real matrix M |25j ; this matrix computes a weight for each term t according to 
the formula Mt. Two well-known orders are lex and grevlex; the former can be 
represented by an identity matrix, and the latter by an upper-triangular matrix 
whose non-zero elements are identical. 

Example 4. Consider the well-known Cyclic-4 system, 

F — (Xl + X 2 + X 3 + X4, X\X 2 + X 2 X 3 + X 3 X4 + X4X1, 
XiX 2 X 3 + X2X3X4 + Z3Z4X1 + X4X1X2, 
XiX 2 X 3 X4 — 1) . 

(1) If we compute a Grobner basis of (F) with respect to lex, we obtain a 
Grobner basis with 6 polynomials made up of 18 distinct terms. 

(2) If we compute a Grobner basis of (F) with respect to grevlex, we obtain a 
Grobner basis with 7 polynomials made up of 24 distinct terms. 

(3) If we compute a Grobner basis of (F) according to the matrix ordering 
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we obtain a Grobner basis with 5 polynomials and 19 distinct terms. 

We can order any finite set of terms using a weight vector in N n , and if necessary, 
we can extend a weight vector to an admissible ordering by adding n — 1 linearly 
independent rows. In the example above, we extended the weight vector (1 3 2 4) 
by adding three more rows. 

The goal of the dynamic algorithm is to discover "good" orderings during the 
Grobner basis computation. Algorithm |2.2| describes a dynamic Buchberger 
algorithm. As with the static algorithm, its basic form contains a number of 
unresolved ambiguities: 

• How shall we select a pair? 

• How do we determine a good ordering? 

• How do we refine the ordering? 

• When should we select a new ordering? 

• Does the algorithm actually terminate? 

2.3. Caboara's implementation. The only implementation of a dynamic algo- 
rithm up to this point was that of 0; it has since been lost. This section reviews 
that work, adapting the original notation to our own, though any differences are 
quite minor. 

2.3.1. How shall we select a pair? Caboara used the sugar strategy, which selects 
(p, q) € P such that the degree of the homogenization of spoly (p, q) is minimal [18 
4 . Pairs were pruned using the Gebauer-Moller algorithm [17] . 

2.3.2. How do we determine a good ordering? Both |20ll9"] suggest using the Hilbert- 
Poincare function to evaluate orderings. Roughly speaking, the Hilbert-Poincare 
function of an ideal I in a ring R, denoted Hr/j (d), indicates: 

• if J is inhomogeneous, the number of elements in R/I of degree no greater 
than d; 
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inputs: F C R 

outputs: G C R and a £ T such that G is a Grobner basis of (F) with respect to 

a 
do: 

(1) Let G = {}, P = {(/,0) : / G P}, a G T 

(2) while P ^ 

(a) Select (p, g) S P and remove it 

(b) Let r be a remainder of spoly (p, g) modulo G 

(c) if r t^O 

(i) Add (g, r) to P for each .g G G 

(ii) Add r to G 

(d) Select r G T 

(e) Let P - P U {(p, q):p,qeG, p^q, lt a (p) # lt T (p)} 

(f) Let cr = r 

(3) return G, cr 

Figure 2.2. Dynamic Buchberger Algorithm 



• if J is homogeneous, the number of elements in R/I of degree d. 
The homogeneous case gives us the useful formula 

H R/I (d) = dim K (R/I) d , 

where the dimension is of the subspace of degree-d terms of the vector space R/I. 
If G is a Grobner basis, then H r //q\ = Pfl/(it(G))i an d it is easy to compute 
the related Hilbert series or Hilbert polynomial for H r //q^ (d) from (It (G)) 
[5J [31 [55]. Many textbooks, such as [52], contain further details. 

In the homogeneous setting, the Hilbert function is an invariant of the ideal 
regardless of the ordering. Thus, we can use it to measure "closeness" of a basis to 
a Grobner basis. Both the static and dynamic algorithms add a polynomial r to 



the basis G if and only if (It (G)) C (It (G U {r})). (See Figure [273]) If we denote 
T = (It (G)) and U = (It (G U {r})), then T C U, so for all d, 

H R/T (d) = dim K (R/T) d > dim K (R/U) d = H R/U (d) . 

If the choice of ordering means that we have two possible values for U, we should 
aim for the ordering whose Hilbert function is smaller in the long run. 

Since G is a Grobner basis only once we complete the algorithm, how can we 
compute the Hilbert function of its ideal? We do not! Instead, we approximate it 
using a tentative Hilbert function H R /n t rQ\\ (d). This usually leads us in the 
right direction, even when the polynomials are inhomogeneous [5]. 

2.3.3. How do we choose the ordering? A potential leading term (PLT) of r G R 
is any term t of r for which there exists an admissible ordering cr such that lto- (r) = 
t. As long as G is finite, there is a finite set of equivalence classes of all monomial 
orderings. We call each equivalence class a cone associated to G, and denote by 
G (cr, G) the cone associated to G that contains a specific ordering a. When g E G, 
we refer to G (cr, t, g) as the cone associated to G that guarantees lt CT (g) = t. 
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FIGURE 2.3. Adding r to G makes Rj (It (G)), and there- 
fore £f(it(G))j smaller. The example here is taken from G = 
{a; 2 + y 2 — 4, xy — 1} with x > y; running the Buchberger algo- 
rithm on the pair ((71,92) gives us It (r) = y 3 , which removes all 
multiples of y 3 + (G) from i?/ (G), thereby decreasing Hn t fQ\\. 

Suppose t £ supp (r) and we want to choose cr such that \t a (r) — t. For each 
u £ supp (r) \ {t}, we want t > u. This means crt > cru, or a (t — u) > 0. A vector 
will suffice to determine cr; and we can find such a vector using the system of linear 
inequalities 

lp (y, t, r) = {y k > 0} n k=1 U J y k (t k - u k ) > | 

(Here, and Uk denote the fcth entry of t and u, respectively.) Since every cone 
G (cr, G) is defined by some system of linear inequalities, this approach successfully 
turns a difficult geometric problem into a well-studied algebraic problem that we 
can solve using techniques from linear programming. 

Proposition 5 (Propositions 1.5, 2.3 of 0). The cone C (cr, F) can be described 
using a union of such systems, one for each f £ F. 

Example 6. In Cyclic-4, the choice of leading terms 

{a;i, 0:10:2, a;ia; 2 2 : 3 I Xia^a^a^} 
can be described by the system of linear inequalities 

{y, > o}ti u {yi - yi > o>t 2 

U {yi - y3 > 0, yi + y 2 - y 3 - y 4 > 0, y 2 - y 4 > 0} 

U {yi - y 4 > 0, y 2 - y 4 > 0, y 3 - y 4 > 0} 

U {yi + y 2 + V3 + 2/4 - > 0} . 

The last inequality comes from the constraint xix 2 X3X 4 > 1, and is useless. This 
illustrates an obvious optimization; if u \ t and the ordering's compatibility 

with division implies that u cannot be a leading term of r; the corresponding linear 
inequality is trivial, and can be ignored. This leads to the first of two criteria to 
eliminate terms that are not potential leading terms. 
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Proposition 7 (The Divisibility Criterion; Proposition 2.5 and Corollary of [9]). 
Let r £ R and t,u€ supp(r). If u divides t properly, then t > T u for every t G T. 
In other words, u is not a potential leading term of r. 

Caboara also proposed a second criterion based on the notion of refining the 
order. If a and r are such that C (r, G)CC (a, G), then we say that r refines er, 
or that r refines the order. Caboara's implementation chooses an ordering r in 



line 2d so that lt r (g) = lt CT (<?) for all g e G\ {r}, so that t refines a. This allows the 
algorithm to discard line |2e] altogether, and makes termination easy to see. While 
this is probably a good idea in general, Golubitsky has found systems where it is 
advantageous to change the ordering, rather than limit oneself to refinement |19| . 

Proposition 8 (The Refining Criterion; Proposition 2.6 of Let G = {gi, . . . , gi\ C 
R and t\, . . . , ti be potential leading terms of g\, . . . , gi, respectively. The system 

{t k > u:ue supp (Sk)}fc=i 

is equivalent to the system 

{t k > u : u a PLT of g k consistent w/tj — It (gj) Vj = 1, . . . , k — l}t. =1 ■ 

Proposition [8] implies that we need only compare tf. with other potential leading 
terms u of gj~ that are consistent with the previous choices; we will call such u, 
compatible leading terms^] From a practical point of view, refining the cone 
in this way is a good idea, as the technique of refining the order allows one to 
warm-start the simplex algorithm from a previous solution using the dual simplex 
algorithm, lessening the overhead of linear programming. It does require some 
record-keeping; namely, retaining and expanding the linear program as we add new 
polynomials. This motivates the definition of 

lp (a, G) = lp (a, {(lt CT (g) , g) : g e G}) := (J lp (y, lt„ (g) , g) . 

gee 

This burden on space is hardly unreasonable, however, as these systems would have 
to be computed even if we allow the order to change. That approach would entail 
a combinatorial explosion. 

Example 9. Let F be the Cyclic-4 system. Suppose that, in the dynamic al- 
gorithm, we add spoly(/i,0) = f\ to G, selecting x\ for the leading term, with 
a = (2, 1, 1, 1). Suppose that next we select spoly (/2, 0) = f2, reduce it modulo G 
to r2 = x\ — 2x2X4 — x% and select x\ as its leading term, with a = (3, 2, 1, 1). We 
have 

lp {a, { (xi,fi) , (xl f 2 ) }) = {y k > 0}£ =1 U { Vl - y k > Q}{ =2 

U {V2 - 2/4 > 0, 2y 2 - 2y 4 > 0} . 

Remark 10. In a practical implementation, it is important to avoid redundant 
constraints; otherwise, the programs quickly grow unwieldy. One way to avoid re- 
dundant constraints is to put them into a canonical form that allows us to avoid 
adding scalar multiples of known constraints. Unfortunately, even this grows un- 
wieldy before too long; one of this paper's main points is to describe a method of 
minimizing the number of required constraints. 



^This notion is essentially Caboara's notion of a potential leading term with respect to F, 
It (F). Our choice of different vocabulary is allows us to emphasize that a "potential" leading term 
for one polynomial is not usually "compatible" with previous choices. 
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FIGURE 2.4. The cones associated to a basis G narrow as we add 
polynomials to the basis. Here, G = {x 2 + y 2 — 4, xy — l}; the 
two possible orderings are r and a, which give lt CT (G) = {x 2 ,xy} 
and lt T (G) = {y 2 ,xy}. Suppose we choose a and compute the 
^-polynomial. This adds r — y 3 +x — Ay to the basis, and the cone 
containing a splits, giving us two choices for It (r). 



2.3.4. When should we refine the ordering? As noted in the introduction, Caboara's 
implementation refines the ordering for a while, then permanently switches the 
refiner off. Making this activation/deactivation mechanism more flexible is the 
major goal of this investigation. 

2.4. The geometric point of view. Here, we provide a visual interpretation of 
how the property G (r, G) C G (cr, G) affects the algorithm. This discussion was 
originally inspired by |20j , but we state it here in terms that are closer to the notion 
of a Grobner fan [23] . 

Any feasible solution to the system of linear inequalities corresponds to a half- 
line in the positive orthant. Thus, the set of all solutions to any given system 
forms an open, convex set that resembles an infinite cone. Adding polynomials to 



G sometimes splits some of the cones. (See Figure 2.4 1. This gives a geometric 



justification for describing Caboara's approach as a narrowing cone algorithm. 



Even though some cones can split when we add new polynomials, not all clones 
must split. In particular, the cone containing the desired ordering need not split, 
especially when the algorithm is nearly complete. For this, reason, it is not neces- 
sary to refine the cone every time a polynomial is added. The methods of the next 
section help detect this situation. 



3. Exploiting the narrowing cone 

The main contribution of this paper is to use the narrowing cone to switch the 
refiner on and off. We propose two techniques to accomplish this: one keeps track 



of cones known to be disjoint (Section 3.1 1; the other keeps track of "boundary 



vectors" that prevent the refiner from leaving the cone (Section 3.2 1 



3.1. Disjoint Cones. The Disjoint Cones Criterion is based on the simple premise 
that if we track inconsistent constraints of linear programs, we can avoid expanding 
the program later. 

3.1.1. Geometric motivation. Let 

• Ce be the cone defined by the selection of t\ , . . . , tt as the leading terms of 
9i, 9i, 
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• C u be the cone defined by the selection of u as the leading term of ge+i, 
and 

• Cg+k be the cone defined by the selection of tg+i, ■ ■ ■ , tg+k as the leading 
terms of g e+1 , g e+k - 

Suppose that a € Cg. If C u n Cg — 0, then it is impossible to refine a to an 
ordering that selects u as the leading term of ge+i, as this would be inconsistent 
with previous choices. 

On the other hand, if r e Cg + k and Cg + k C Cg, then it is possible to refine a to 
an ordering r. Let 

• C v be the cone defined by the selection of v as the leading term for gg+k+i- 

If C v C C u , then (7„ n Cg+k Q C U C\ Cg — 0, so we cannot refine a to any ordering 
that selects v as the leading term of gg+k+i- Is there some way to ensure that the 
algorithm does not waste time trying to solve such systems of linear inequalities? 
Yes! If we record the cone C u , we can check whether C v C C u , rather than going 
to the expense of building a linear program to check whether C v n Cg+k 0- 

3.1.2. Algebraic implementation. Determining accurately that C v C C„ is not so 
easy; we will content ourselves with determining when one set of inconsistent linear 
constraints is a subset of another set. 

Theorem 11 (Disjoint Cones Criterion). Let L u , L v , and L a be sets of linear 
constraints. If L a is inconsistent with L u and L u C L v , then L v is also inconsistent 
with L a . 

Proof. Assume L a is inconsistent with L u and L u C L v . The first hypothesis 
implies that L a n L u = 0. The second implies that L v has at least as many 
constraints as L u , so that the feasible regions C u and C v , corresponding to L u and 
L v , respectively, satisfy the relation C u D C v . Putting it all together, C v n C a C 

C u nC a = <D. □ 

In our situation, L u and L v correspond to different choices of leading terms for 
polynomials in the set, while L a = lp (a, G). We will want to consider both the case 
where L v is the set of new constraints, and the case where L v is some extension of 
L a . Rather than discard the inconsistent linear program L u , we will retain it and 
test it against subsequent sets of constraints, avoiding pointless invocations of the 
simplex algorithm. 

We implement the geometric idea using a global variable, rejects. This is a set of 
sets; whenever L a is known to be consistent, but the simplex algorithm finds L a L)L u 
inconsistent, we add L u to rejects. Subsequently, while creating the constraints in 
an extension L v of L a , we check whether L u is contained in these constraints; if so, 
we reject L v out of hand, without incurring the burden of the simplex algorithm. 

Again, we are not checking whether the cones are disjoint, only whether the 
linear constraints are a subset. With appropriate data structures, the complexity 
of determining subset membership is relatively small; with hashed sets, for exam- 
ple, the worst-case time complexity would be the cost of the hash function plus 
O (\rejects\). 

3.2. Boundary Vectors. Unlike the Disjoin Cones Criterion, the Boundary Vec- 
tors Criterion can prevent the construction of any constraints. 
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3.2.1. Geometric motivation. Let C (G,cr) be a cone associated with G, containing 
the ordering a. 

Definition 12. The closure of C (G, cr) is the feasible region obtained by including 
the boundaries of the solution of lp (cr,G). We say that oj G K™ is a boundary 
vector of G (G, cr) if 

• n = 2 and w lies on the closure of G (G, <j), or 

• ?i > 2 and oj lies on the intersection of two non-redundant hyperplanes that 
define G(G,ct). 

We denote the set of all boundary vectors of G (G, cr) by £Ig,v If G and cr are 
understood from context, we simply write f2. 

Example 13. Suppose a = (18, 5, 7) and G (G, cr) is defined by 

2yi - vi > 
-2/i + 4y 2 > 
2/i + 2/2 - 3y 3 > 
-2/2 +2/3 > 

We can find D, = {(68, 17, 17) , (204, 51, 85) , (34, 17, 17)} by taking intersections 
of the hyperplanes that define the cone. However, not all such intersections are 
boundary vectors; two intersections of the hyperplanes do not border the feasible 
region are (17,34,17) and (17,34,34). 

Boundary vectors possess the very desirable property of capturing all possible 
refinements of the term ordering. 

Theorem 14 (Boundary Vectors Criterion). Let r G R and f2 the set of boundary 
vectors of C (cr, G) . Write t = lt CT (r) . If there exists r G G (cr, G) such that lt T (r) = 
u =/= t — that is, there exists r that refines the order — then there exists oj G Q 
such that ui (u — t) > 0. 

Two observations are in order before we prove Theorem |14| First, the converse 
of Theorem 14 is not true in general. For example, let oj — (2, 1), r = x 2 + x + y, 



v = x 2 . u — x, and t — y. Even though u> (u — t) > 0, the fact that u \ v implies 
that no admissible ordering r chooses lt T (r) = u. 

Nevertheless, the theorem does imply a useful corollary: 

Corollary 15. // we know the set Q of boundary vectors for C (cr, G), then we can 
discard any term u such that Li (t — u) > for all w £ O. That is, u is not a 
compatible leading term. 



We turn now to the proof of Theorem [14] Figure |3.1| illustrates the intuition; 
here, r and a select different leading terms, and oj lies on the other side of r from 
cr. By linearity, oj gives greater weight to lt T (r) than to lt CT (r). 



[of Theorem 14 Suppose that there exists r € G (cr, G) so that lt T (r) = u. By 
definition, r (u — t) > 0. If r = oj for some u G O, then we are done. Oth- 
erwise, consider the linear program defined by maximizing the objective function 
(y%, ... , y n ) ■ (u — t) subject to the closure of L — lp (cr, G) U {y% + ■ ■ ■ + y n = 1}. 
This is a convex set; the well-known Corner Point Theorem implies that a maxi- 
mum occurs at a corner point of the set [TH [28] : by definition, this solution is a 
boundary vector of G (cr, G). Let oj G G (cr, G) be a boundary vector where L takes 
its maximum value; then oj (u — t) > r (u — t) > 0. 
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FIGURE 3.1. If lt T (r) 7^ lt CT (r), convexity and linearity imply that 
we can find a boundary vector uj such that uj would give lt r (r) more 
weight than lt^ (r) . 



Example 16. We illustrate the linear program in the proof using Example [13] 
Rewrite the strict inequalities to allow equality, 



2/1 



2yi - y2 

-V2 + M)2 

2/2 - 3t/ 3 

-2/2 + J/3 



> o 

> o 

> o 

> 



then add the constraint 



2/1 



•3/2 + 3/3 = 1. 

We obtain a two-dimensional cross-section of C (er, G) . This cross-section has corner 
points 

'2 3 3 \ /3 3 1\ /l 1 r 
3' 18' 18/ ' V 5 ' 20' 4 / ' \2'4'4 



which scale to the boundary vectors that we found in Example 13 



Computing Q can be impractical, as it is possibly exponential in size. We ap- 
proximate it instead by computing corner points that correspond to the maximum 
and minimum of each variable on a cross section of the cone with a hyperplane, 
giving us at most 2n points. Figure [3~2| illustrates the idea. 



Example 17. Continuing Example |16| maximizing yi, y 2 , and 2/3 gives us the 
boundary vectors listed in that example. 

Approximating 51 has its own disadvantage; inasmuch as some orderings are 
excluded, we risk missing some refinements that could produce systems that we 
want. We will see that this is not a serious drawback in practice. 



3.2.2. Minimizing the number of constraints. While the optimization of Section [3T2 
reduces the number of linear programs computed, their size can remain formidable 
all the same. This is because we are still adding constraints for every monomial 
that passes the Divisibility Criterion. Is there some way to use boundary vectors 
to minimize the number of constraints in the program? 
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FIGURE 3.2. We approximate fl by computing boundary vectors 
corresponding to points that maximum and minimize the value of 
an objective function. The cone at the left has seven boundary 
vectors, as we see in the cross section on the right. In this case, 
four vectors maximize and minimize the variables; they define a 
"sub-cone" whose cross-section corresponds to the dashed line. 



/<* * \ / V 




(a) (b) 

Figure 3.3. Although a and /i lie within the same cone, the choice 
of border vectors in (a) could mean the algorithm at first does not 
add a constraint to guarantee lt CT (gi) > lt M (gi). If the narrowed 
cone splits later on, as in (b), and the algorithm moves into a sub- 
cone that does not contain a, a subsequent choice of /i is possible, 
causing a change in the leading terms. We can try adding the 
previously-overlooked constraint, and continue if we find a feasible 
solution. The unlabeled dot could represent such a compromise 
ordering. 



We will attempt to add only those constraints that correspond to terms that the 
boundary vectors identify as compatible leading terms. As we are not computing 
all the boundary vectors, the alert reader may wonder whether this is safe. 

Example 18. Suppose 

• /i E C (r, {gi, . . .,ge+k}) C C (<r,{gi, . . .,&}), 

• t = Iter (gt), and 

• u = lt p (gi). 

Suppose further that, when gi is added to the basis, the algorithm selects a for the 
ordering. 

For some choices of boundary vectors, the algorithm might not notice that 
(J, € C (a, {g\, . . . ,<%}), as in Figure [3~3^ a). Thus, it would not add the constraint 
(y%, . . . ,y n ) • (t — u) > 0. A later choice of boundary vectors does recognize that 
fi G C (r, {gi, . . . , gi + k}), and even selects fi as the ordering. We can see this in 
Figure [3~3j [b) . In this case, the leading term of gi changes from t to u; the ordering 
has been changed, not refined! 
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Figure 3.4. Estimation of the corner points of C (<r,G) by a 
perturbation of the inequalities. Each diagram shows a two- 
dimensional cross-section of a three-dimensional cone. The dashed 
lines on the connect actual boundary vectors. The diagram on the 
right connects boundary vectors of a perturbation of the system of 
linear inequalities. The perturbation is designed to give us interior 
points of the cone that are close to corner points. 

Since the Grobner basis property depends on the value of the leading terms, 
this endangers the algorithm's correctness. Fortunately, it is easy to detect this 
situation; the algorithm simply monitors the leading terms. 

It is not so easy to remedy the situation once we detect it. One approach 
is to add the relevant critical pairs to P, and erase any record of critical pairs 
computed with g£, or discarded because of it. Any practical implementation of the 
dynamic algorithm would use criteria to discard useless critical pairs, such as those 
of Buchberger, but if the polynomials' leading terms have changed, the criteria no 
longer apply. How do we recover those pairs? 

A second approach avoids these quandaries: simply add the missing constraint 
(y%, . . . ,y n ) • (t — u) > to the system. This allows us to determine whether we 
were simply unlucky enough to choose a from a region of C (er, G) that lies outside 
C (/i, G), or whether really there is no way to choose both lt CT (ge) and lt M (ge+k) 
simultaneously. In the former case, the linear program will become infeasible, and 
we reject the choice of ^; in the latter, we will be able to find t e C (er, G)nC (/x, G). 

3.2.3. Implementation. Since the constraints of lp (<r, G) consist of integer polyno- 
mials, it is possible to find integer solutions for the boundary vectors; one simply 
rescales rational solutions once they are found. However, working with exact arith- 
metic can be quite slow, techniques of integer programming are very slow, and for 
our purposes, floating-point approximations are quite suitable. Besides, most linear 
solvers work with floating point numbers. 

On the other hand, using floating point introduces a problem when comparing 
terms. The computer will sometimes infer uj ■ (u — t) > even though the exact 
representation would have uj ■ (u — t) = 0. We can get around this by modifying 
the constraints of the linear program to uj ■ (u — t) > e for some sufficiently large 
e > 0. As we see in Figure |3~4] the polygon no longer connects extrema, but points 
that approximate them. We might actually reject some potential leading terms u 
on this account, but on the other hand, we never waste time with terms that are 
incompatible with the current ordering. 

This modified linear program is in fact useful for computing feasible points to 
the original linear program as well. 
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Theorem 19. Let e > and J a finite subset o/N. The system of linear inequalities 
lp (a, G) = |a« • (yi, . . . , j/„) > o} U > 0}" 



is feasible if and only if the linear program 

{a«.( tfl ,...,y n )>e} U { Vj > e}^ 

is also feasible. 

A solution to the linear program obviously solves the system of linear inequalities. 
Thus, suppose the system of linear inequalities has a solution a . Let r = ba, where 
b £ N is chosen large enough that tx- > e for k = 1, . . . , n. For each j £ J, define 



fc=i 



then choose Cj > 1 such that Cj7j > e. Put = max{cj} jg 7 . Let w = dr. We have 
Wfc > Tfc > e for each k = 1, . . . , n, and 



fe=i fc=i 

for each j S J. We have shown that uj is a solution to the linear program, which 
means that the linear program is also feasible. 

Based on Theorem [19] we take the following approach: 

(1) Replace each constraint eft' ■ y > of lp (a, G) with afj> ■ y > e, add 
constraints y^ > e for fc = 1, . . . ,n, and take as the objective function the 
minimization of ^ t/fc. We denote this new linear program as mlp (a, G). 

(2) Solve mlp (a,G). This gives us a vector t that can serve as a weighted 
ordering for the terms already computed. 

(3) Identify some d £ R such that mlp (<r, G) intersects the hyperplane Vk = 
d, giving us a cross-section X of the feasible region. This is trivial once we 
have a solution r to mlp {a, G), since we can put d = 1 + T k- 

(4) Compute an approximation to Q by maximizing and minimizing each y^ 
on K. 



Algorithm compute_ boundary _ vectors ( Figure 3.5) gives pseudocode to do this; it 
generates a set of boundary vectors ^ that approximates the set f2 of boundary 
vectors of C (a, G). 

Once we have an approximation VP to the boundary vectors f2, we use it to elim- 
inate terms that cannot serve as leading terms within the current cone. Algorithm 



identify _ cits _using_boundary_vectors (Figure 3.6 I, accomplishes this by looking 
for u £ supp (r) \ {lto- (r)} and ip £ such that if) (u — t) > 0. If it finds one, then 
u is returned as a compatible leading term. 

In Section |3.2.2| we pointed out that creating constraints only for the terms 
identified as potential leading terms by the use of boundary vectors can lead to 
an inconsistency with previously-chosen terms. For this reason, we not only try to 



solve the linear program, but invokes algorithm monitor _lts (Figure 3.7 1 to verify 
that previously-determined leading terms remain invariant. If some leading terms 
would change, the algorithm obtains a compromise ordering whenever one exists. 

Theorem 20. Algorithm monitor_lts of Figure \3. 7| terminates correctly. 
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algorithm compute_ boundary _ vectors 
inputs: 

• L — mlp (a, G) 

• t eW l that solves L 

outputs: $C1°, a set of vectors that approximates the boundary vectors of 

mlp (cr, G) 
do: 

(1) let * = {} 

(2) let d = 1 + n + • • • + r„ 

(3) let L = L U { yi + ■ ■ ■ + y n = d} 

(4) for fce{l,...,n} 

(a) add to the solution w of L that maximizes ujk 

(b) add to ^ the solution w of I that minimizes 

(5) return ^ 

Figure 3.5. Algorithm to compute approximate boundary vec- 
tors to C (G, cr) 



algorithm identify _ clts_ using_ boundary _ vectors 
inputs: 

• cr e T, the current term ordering 

• t = \t a (r), where r e R 

• U — supp (r) \ {t} 

• ^ C E™, approximations to the boundary vectors of C (G,a) 
outputs: V C [/, where ti€Fiffti = tor«eC/ and V (v — t) > for some 

V> e * 

do: 

(1) let V = {t} 

(2) for u € C/ 

(a) if -0 (u — t) > for some ip £ ^ 
(i) add u to V 

(3) return V 

FIGURE 3.6. Eliminating terms using approximation to boundary vectors 



Proof. Termination is evident from the fact that G is a finite list of polynomials, 
so the while loop can add only finitely many constraints to L. Correctness follows 
from the fact that the algorithm adds constraints to mlp (t, G) if and only if they 
correct changes of the leading term. Thus, it returns (True, /x) if and only if it 
is possible to build a linear program L whose solution [i lies in the non-empty set 
C(T,G)nC(or,G). * □ 
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algorithm monitor _ Its 
inputs 

• G, the working basis 

• a E T, the old ordering 

• r E T, a new ordering 

• L = mlp (t, G) 

outputs 

• (True, /i) if there exists fj, E L refining G (G,<r) and lt M (giast) = lt T (<7iast)> 
where <?i ast is the newest element of G 

• False otherwise 

do 

(1) let fj, = t 

(2) while there exists g E G such that lt M (g) ^ lt CT (g) 

(a) for each g E G whose leading term changes 

(i) let t = lt CT (», u = \t ll {g) 

(ii) letL = LU{y(t-u)} 

(b) if L is infeasible return False 
else let be the solution to L 

(3) return (True, /x) 

FIGURE 3.7. Ensuring the terms remain invariant 



4. Experimental results 

The current study implementation, written in Cython for the Sage computer 
algebra system [57], is available at 

www . math . usm . edu/perry/Resear ch/dynamic_gb . pyx 
It is structured primarily by the following functions: 

dynamic_gb is the control program, which invokes the usual functions for a Buch- 
berger algorithm (creation, pruning, and selection of critical pairs using the 
Gebauer-Moller algorithm and the sugar strategy, as well as computation and 
reduction of of ^-polynomials) , as well as the following functions necessary for 
a dynamic algorithm that uses the criteria of Disjoint Cones and Boundary 
Vectors: 

choose_an_ordering, which refines the ordering according to the Hilbert 
function heuristic, and invokes: 

possible_lts, which applies the Boundary Vectors and Divisibility cri- 
teria; 

feasible, which tries to extend the current linear program with con- 
straints corresponding to the preferred leading term; it also applies 
the Disjoint Cones criterion, and invokes 

monitor_lts, which verifies that an ordering computed by feasible 
preserves the previous choices of leading terms; 
boundary_vectors, which computes an approximation ^ to the boundary 
vectors ft. 

The dynamic_gb function accepts the following options: 
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• static: boolean, True computes by the static method, while False (the 
default) computes by the dynamic method; 

• strategy: one of 'sugar', 'normal' (the default), or 'mindeg'; 

• weighted_sugar: boolean, True computes sugar according to ordering, 
while False (the default) computes sugar according to exponents; 

• use_boundary_vectors: boolean, default is True; 

• use_disjoint_cones: boolean, default is True. 

At the present time, we are interested in structural data rather than timings. To 
that end, experimental data must establish that the methods proposed satisfy the 
stated aim of reducing the size and number of linear programs constructed; in other 
words, the algorithm invokes the refiner only when it has high certainty that it is 
needed. Evidence for this would appear as the number of linear programs it does 
not construct, the number that it does, and the ratio of one to the other. We should 
also observe a relatively low number of failed linear programs. 

Table [T] summarizes the performance of this implementation on several bench- 
marks. Its columns indicate: 

• the name of a polynomial system tested; 

• the number of linear programs (i.e., potential leading terms) 

— rejected using approximate boundary vectors, 

— rejected using disjoint cones, 

— solved while using the two new criteria, 

— solved while not using the new criteria, and 

— failed; 

• the ratio of the linear programs solved using the new criteria to the number 
solved without them; and 

• the number of constraints in the final program, both using the new criteria, 
and not using them. 

Both to emphasize that the algorithm really is dynamic, and to compare with 
Caboara's original results, Table [T] also compares: 

• the size of the Grobner basis generated by the dynamic algorithm, in terms 
of 

— the number of polynomials computed, and 

— the number of terms appearing in the polynomials of the basis; 

• the size of the Grobner basis generated by Singular's std() function with 
the grevlex ordering, in the same terms. 

The systems "Caboara i" correspond to the example systems "Es i" from [9], some 
of which came from other sources; we do not repeat the details here. We verified 
by brute force that the final result was a Grobner basis. 

The reader readily sees that the optimizations introduced in this paper accom- 
plish the stated goals. By itself, the method of boundary vectors eliminates the 
majority of incompatible leading terms. In the case of dense polynomial systems, 
it eliminates the vast majority. This means that far, far fewer linear programs are 
constructed, and those that are constructed have far fewer constraints than they 
would otherwise. While the number of inconsistent linear programs the algorithm 
attempted to solve (the "failed" column) may seem high in proportion to the number 
of programs it did solve ("solved") , this pales in comparison to how many it would 
have attempted without the use of boundary vectors; a glance at the "div only" 
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*This system was terminated when using only the Divisibility Criterion, as the 
linear program had acquired more than 1000 constraints. 
Table 1 . Dynamic algorithm with sugar strategy, applying Divis- 
ibility Criterion after boundary vectors. Data for static algorithm 
included for comparison. 
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column, which consists of monomials that were eliminated only by the Divisibility 
Criterion, should allay any such concerns. 

Remark 21. In two cases, the new criteria led the algorithm to compute more linear 
programs than if it had used the Divisibility Criterion alone. There are two reasons 
for this. 

• One of the input systems (Caboara 2) consists exclusively of inhomoge- 
neous binomials with many divisible terms. This setting favors the Divisi- 
bility Criterion. The other system (Caboara 5) also contains many divisible 
terms. 

• For both systems, the sample set of boundary vectors wrongly eliminates 
compatible monomials that should be kept. While this suggests that the 
current strategy of selecting a sample set of boundary vectors leaves much 
to be desired, the dynamic algorithm computes a smaller basis than the 
static, and with fewer iS-polynomials, even in these cases. 

It should not startle the reader that the dynamic algorithm performs poorly on 
the homogeneous Katsura-n systems, as Caboara had already reported this. All 
the same, boundary vectors and disjoint cones minimize the cost of the dynamic 
approach. 

Many of our Grobner bases have different sizes from those originally reported by 
Caboara. There are several likely causes: 

• The original report had typographical errors. We have verified and cor- 
rected this in some cases, but some results continue to differ. 

• The static ordering used here may order the variables differently from [9], 
which did not documented this detail. For example, Caboara reports a basis 
of only 318 polynomials for Caboara 2 when using the static algorithm with 
grevlex, but Sage (using Singular) finds 553. 

• Several choices of leading term can have the same tentative Hilbert function, 
but one of the choices is in fact better than the others in the long run. 
The original implementation may have chosen differently from this one. In 
particular, this implementation applies the Hilbert heuristic to select an 
initial ordering only on the first polynomial, while [9 applied it to more 
than one. 

We conclude this section with a word on complexity. While the worst case time 
complexity of the simplex algorithm is exponential |21| . on average it outperforms 
algorithms with polynomial time complexity. The computation of boundary vec- 
tors increases the number of invocations of simplex, but the large reduction in the 
number of monomials considered more than compensates for this. The space re- 
quirements are negligible, as we need only 2n boundary vectors at any one time, 
and the reduction in the number and size of the linear programs means the algo- 
rithm needs to remember only a very few disjoint cones. Considering how rarely 
the disjoint cones are useful, it might be worthwhile not to implement them at all, 
but we have not observed them to be a heavy burden at the current time. 

5. Conclusion, future work 

We set out to reduce the size and number of linear programs used by a dynamic 
algorithm to compute a Grobner basis. Geometrical intuition led us to two methods 
that work effectively and efficiently. While the effect with the systems tested by 
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Caboara was only moderate, and in some cases counterproductive, the story was 
different with the dense benchmark systems. In these cases, the number of refine- 
ments approached insignificance, and we continued to achieve good results. The 
final Grobner basis was always of a size significantly smaller than grevlex; only a 
few refinements were required for the algorithm to be very effective. 

A significant restraint imposed by many computer algebra systems is that a term 
ordering be defined by integers. As Sage is among these, this has required us to 
solve not merely linear programs, but pure integer programs. Integer programming 
is much more intensive than linear programming, and its effect was a real drag on 
some systems. There is no theoretical need for this; we plan to look for ways to 
eliminate this requirement, or at least mitigate it. 

An obvious next step is to study various ambiguities in this approach. This 
includes traditional questions, such as the effect of the selection strategy of critical 
pairs, and also newer questions, such as generating the sample set of boundary 
vectors. We followed Caboara's approach and used the sugar strategy. We expect 
the normal strategy [5] to be useful only rarely, while signature-based strategies |15| 
[13], which eliminate useless pairs using information contained in the leading terms 
of a module representation, could combine with the dynamic algorithm to expand 
significantly the frontiers of the computation of Grobner bases, and we are working 
towards such an approach. 

Another optimization would be to use sparse matrix techniques to improve the 
efficiency of polynomial reduction in a Grobner basis algorithm, a la F4. Combining 
this with a dynamic algorithm is comparable to allowing some column swaps in the 
Macaulay matrix, something that is ordinarily impossible in the middle of a Grobner 
basis computation. 

The authors would like to thank Nathann Cohen for some stimulating conversa- 
tions, and his assistance with Sage's linear programming facilities. 

Appendix 

Here we give additional information on one run of the inhomogencous Cyclic-6 
system. Computing a Grobner basis for this ideal with the standard sugar strategy 
requires more S'-polynomials than doing so for the ideal of the homogenized system; 
in our implementation, the number of S'-polynomials is roughly equal to computing 
the basis using a static approach. (In general, the dynamic approach takes fewer 
S-polynomials, and the weighted sugar strategy is more efficient on this ideal.) 
Information on timings was obtained using Sage's profiler: 

• for the static run, we used 

7,prun B = dynamic_gb(F,static=True,strategy='sugar'); 

• for the dynamic run, we used 

7,prun B = dynamic_gb(F,static=False,strategy='sugar'). 

We obtained structural data by keeping statistics during a sample run of the 
program. 

While Cython translates Python instructions to C code, then compiles the result, 
the resulting binary code relies on both the Python runtime and Python data 
structures; it just works with them from compiled code. Timings will reflect this; 
we include them merely to reassure the reader that this approach shows promise. 
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Our implementation of the dynamic algorithm uses 377 5-polynomials to com- 
pute a Grobner basis of 20 polynomials, with 250 reductions to zero. At termination. 
rejects contained 43 sets of constraints, which were never used to eliminate linear 
programs. As for timings, one execution of the current implementation took 25.45 
seconds. The profiler identified the following functions as being the most expensive. 

• Roughly half the time, 11.193 seconds, was spent reducing polynomials. 
(All timings of functions are cumulative: that is, time spent in this func- 
tion and any other functions it invokes.) This is due to our having to 
implement manually a routine that would reduce polynomials and compute 
the resulting sugar. A lot of interpreted code is involved in this routine. 

• Another third of the time, 8.581 seconds, was spent updating the critical 
pairs using the Gebauer-Moller update algorithm. Much of this time (3.248 
seconds) was spent computing the 1cm of critical pairs. 

The majority of time (5/6) was spent on functions that are standard in the Buch- 
berger algorithm! Ordinarily, we would not expect an implementation to spend 
that much time reducing polynomials and updating critical pairs; we are observing 
a penalty from the Python runtime and data structures. The size of this penalty 
naturally varies throughout the functions, but this gives the reader an idea of how 
large it can be. 

Other expensive functions include: 

• About 4.6 seconds were spent in choose_an_ordering. Much of this time 
(2.4 seconds) was spent sorting compatible leading terms according to the 
Hilbert heuristic. 

• The 238 invocations of Sage's hilbert_series and hilbert_polynomial 

took 1.5 seconds and 0.8 seconds, respectively. As Sage uses SINGULAR 
for this, some of this penalty is due to compiled code, but Singular's 
implementation of these functions is not state-of-the-art; techniques derived 
from [5] and |26| would compute the Hilbert polynomial incrementally and 
quickly. 

• The 57 invocations of feasible consumed roughly 1.2 seconds. This in- 
cludes solving both continuous and pure integer programs, and suffers the 
same penalty from interpreted code as other functions. 

It is worth pointing out what is not present on this list: 

• Roughly one quarter of a second was spent in monitor_lts. 

• The boundary_vectors function took roughly one and a half tenths of a 
second (.016). 

• Less than one tenth of a second was spent applying the Boundary Vectors 
Criterion in possible_lts. 

If we remove the use of the Boundary Vector and Disjoint Cones criteria, the relative 
efficiency of feasible evaporates; the following invocation illustrates this vividly: 

7«prun B = dynamic_gb(F , strategy= ' sugar ' , 
use_boundary_vectors=False , 
use_disjoint_cones=False) 

Not only is the Divisibility Criterion unable to stop us from solving 2,820 linear 
programs, but the programs themselves grow to a size of 247 constraints. At 41 
seconds, feasible takes nearly twice as long as entire the computation when us- 
ing the new criteria! This inability to eliminate useless terms effect cascades; 
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hilbert_polynomial and hilbert_series are invoked 3,506 times, taking 12.6 
and 28.6 seconds, respectively; choose_an_ordering jumps to 92.7 seconds, with 
44 seconds wasted in applying the heuristic. By contrast, the timings for reduction 
of polynomials and sorting of critical pairs grow much, much less, to 20 seconds and 
12.9 seconds, respectively. (The increase in reduction corresponds to an increase 
in the number of ^-polynomials, probably due to a different choice of leading term 
at some point - recall that the approximation of boundary vectors excludes some 
choices.) The optimizations presented here really do remove one of the major bot- 
tlenecks of this method. 

We conclude by noting that when we use the Disjoint Cones criterion alone, 
1,297 invalid monomials are eliminated, but the number of constraints in the final 
linear program increases to 250; this compares to 4,080 monomials that Boundary 
Vectors alone eliminates, with 56 constraints in the final linear program. 
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